Astronomy & Astrophysics manuscript no. 1415 1 
May 6, 2010 



©ESO2010 



Spatially resolved XMM-Newton analysis and a model of the 
nonthermal emission of MSH15-52 

KM. Schock 1 , 1. Biisching 2 , O.C. de Jager 2 , P. Eger 1 , and M.J. Vorster 2 



1 Erlangen Centre for Astroparticle Physics, Erwin-Rommel-StraBe 1, 91058 Erlangen, Germany 

2 Unit for Space Physics, North-West University, Potchefstroom, 2520, South Africa 



Preprint online version: May 6, 2010 



ABSTRACT 



o 



in 

< 

o 

Or 
I ' 

o- 



> 

in 

00 

q 

in 

o 
o 



X 



We present an X-ray analysis and a model of the nonthermal emission of the pulsar wind nebula (PWN) MSH 15-52. We analyzed 
XMM-Newton data to obtain the spatially resolved spectral parameters around the pulsar PSRB1509-58. A steepening of the fitted 
power-law spectra and decrease in the surface brightness is observed with increasing distance from the pulsar. In the second part of 
this paper, we introduce a model for the nonthermal emission, based on assuming the ideal magnetohydrodynamic limit. This model 
is used to constrain the parameters of the termination shock and the bulk velocity of the leptons in the PWN. Our model is able to 
reproduce the spatial variation of the X-ray spectra. The parameter ranges that we found agree well with the parameter estimates found 
by other authors with different approaches. In the last part of this paper, we calculate the inverse Compton emission from our model 
and compare it to the emission detected with the H.E.S.S. telescope system. Our model is able to reproduce the flux level observed 
with H.E.S.S., but not the spectral shape of the observed TeV y-ray emission. 
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1. Introduction 

MSH 15-52 (G320.4-1.2) is a complex supernov a remnant 

iSNR) that was discovered in radio wavelengths by iMills et all 
1961). Its radio appearance is dominated by two spots, a 
brighter one to the northwest and a fainter one to the south- 
east. The radio emission to the northwest coincides spatially with 
the HII region RCW89. Situated inside the SNR is the ener- 
getic 150 ms pulsar PS RB1509-58, which wa s discovered by 
the Einstein satellite dSeward & Harndenlll982l) . It has a spin- 
down luminosi ty of £ - 1.8 X 10 3 7 erg/s and a characteristic age 
of T = 1 .7 kyr (Kas pi et al1[l9 94). This makes the pulsar one of 
the youngest and most energetic known. Gaens ler et al.l (1999) 
concluded by a comparison of the radio and the X-ray emission 
that MSH 15-52, PSRB1509-58 and RCW89 are associated 
objects and are located at a distance of (5.2 ± 1 .4) kpc. 

X-ray observations of the SNR have sh own a pulsar 
wind nebula (PWN) powered by PSRB 1509-58 dTrussoni et all 
1996; Tamura et al. 1996). Observations with the Chandra satel- 
lite have revealed two outflow jets in the southeast and northwest 
directions, the latter t erminating in the optical nebula RCW 89 
dGaensleretal .1120021 henceforth referred to as G02). G02 de- 
rived a power-law photon index of 2.05 ± 0.04 and an absorp- 
tion column density of (9.5 ± 0.3) xlO 21 cm -2 for the diffuse 
PWN in the energy band of 0.5-10keV. A detailed study of the 
innerm ost region of MSH 15-52 was performed by I Yatsu et al.l 
(l2009l henceforth referred to as Y09) using an extended data 
set of Chandra observations. Their analysis revealed hints for a 
ring-like feature around the pulsar, which might correspond to a 
wind termination shock. The hard X-ray spectrum of the PWN 
was observed with the BeppoSAX and the INTEGRAL satel- 
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Fig. 1. XMM-Newton count map of MSH 15-52 showing the six 
extraction regions (numbered and marked by white lines) and the 
exclusion region (black dashed line) used for the spectral anal- 
ysis. The six extraction regions are centered on the position of 
the pulsar PSRB 1509-58. The exclusion region was chosen to 
encompass the thermal emission from the H II region RCW 89. 



lites. The PWN is clearly seen in the off-pulse component of the 
emission and the morphology corresponds nicely to the measure- 
ments at lower energies. In the energy band of 20-200 keV, the 
off-pulse emission from the PWN is fitted best by a power law 
with an index of 2.1 dMineo et aLfe OOl: Forot et al. 2006). 

In the very high-energy (VHE; 100 GeV-100 TeV) y-ray do- 
main, the PWN was observed by the High Energy Stereoscopic 
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Table 1. Details of the XMM-Newton EPIC PN observations on 
MSH 15-52. 



Observation 
ID 



Exposures (ks) 
performed (1> net 



(2) 



0207050201 
0302730201 
0302730301 



23.135 
16.130 

8.235 



5.9 
3.6 
2.0 



(1) Exposure time without background screening 

(2) Net exposure time after background screening 

Table 2. Extraction regions for the XMM-Newton spectral anal- 
ysis of the source. 



Ring 


Radius ( 


arcsec) 


No. ofObs. 


No. 


inner 


outer 




1 


30 


57 


2 


2 


57 


84 


2 


3 


84 


138 


3 


4 


138 


192 


3 


5 


192 


246 


1 


6 


246 


300 


2 



blocking filters. The observations used in our analysis are listed 
in Table[T] 

For the analysis of the X-ray data we used the XMM-Newton 
Science Analysis System (SAS) version 8.0.0 supported by tools 
from the FTOOLS package. For the spectral modeling, version 
12.5.0 of the XSPEC software was used (Arnaud 1996). To 
screen the data from periods of high background-flaring activ- 
ity, we used the 7 to 15 keV lightcurve provided by the standard 
SAS analysis chain extracted from the full FoV. Since a good 
understanding of the background is crucial for the analysis of 
extended sources, we applied a conservative background thresh- 
old of ten background counts per second for the definition of the 
good time intervals (GTI). Furthermore, we only analyzed the 
data of the PN camera since it is more sensitive than the MOS 
cameras. We selected good (FLAG= 0) single and double events 
(PATTERN<= 4). The innermost region used for the spectral 
analysis starts at a radius of 30 arcsec from the bright pulsar and 
therefore, pile-up is not an issue. All of the observations were 
affected by long periods of background flaring or full scientific 
buffer of the PN camera. This leads to rather short net exposures 
(see TableQ}. However, the statistics are still sufficient to obtain 
spectra from all of the extraction regions as defined in the next 
section. 



System (H.E.S.S.). The source is clearly extended beyond the 
point spread function (PSF) and shows a morphology compara- 
ble to what is seen in X-rays, extending northwest and southeast 
of the pulsar. The observed y-ray emission is well-fitted by a 
power l aw with a photon index of 2.27 up to a photon energy of 
40 TeV dAharonian et al.ll2005h . 

In the first part of this paper, we present an analysis of the 
XMM-Newton data of the PWN MSH 15-52. Its large effective 
area makes XMM-Newton ideally suited for the spectral analysis 
of extended sources. Our analysis thus provides a good measure- 
ment of the large-scale characteristics of MSH 15-52, compared 
to earlier high-resolution measurements of the inner region with 
the Chandra satellite. Following the XMM-Newton analysis, we 
introduce a model which is based on the assumption of the ideal 
magnetohydrodynamic limit to describe the observed emission 
in the X-ray domain. Fitting the spatially resolved X-ray emis- 
sion with the model, we found optimum parameters for several 
scenarios that have been discussed by other authors. This allows 
us to constrain physical quantities of the PWN termination shock 
and the velocity profile of the PWN. In the last part, we also ap- 
ply this model to make predictions for the VHE y-ray emission 
as observed by the H.E.S.S. Cherenkov telescope array. 

2. XMM-Newton Observations 

The region around the SNR MSH 15-52 has b een observed six 
times with XMM-Newton dJansen et all 1200 ll) with the EPIC- 
MOS (lTurner et alJ l200lt) and EPIC-PN (tStriider et al.1 l200l 
cameras. These pointings were either centered north or south- 
east of the pulsar PSR B 1509-58. During one of these pointings 
(Observation ID: 0312590101), all three cameras were operated 
in timing mode. Therefore, this observation will not be consid- 
ered in the present paper. In order to study the morphology- 
dependent spectral characteristics of the extended emission, we 
require the whole area of our interest to be within the field of 
view (FoV). Only the three observations pointing towards the 
northern region matched this criterion. In each of them the de- 
tectors were operated in full-frame mode with medium optical 



3. X-ray Spectra 

We extracted spectra from annular regions centered on the pulsar 
position (R. A.: 15:13:55,Dec: -59:08:08.8) to study the changes 
of the spectral properties. The regions can be seen on the XMM- 
Newton sky map (Fig.[TJ. We use full annuli in contrast to wedge 
shapes, to extract the integral flux for each radial distance. This is 
especially important for the modeling of the flux of the source, 
which is presented in Section [4] The parameters for the rings 
are listed in Table [2] The numbering scheme starts with ring 1 , 
which is closest to the pulsar, and goes out to ring 6, for which 
the outer radius lies at a distance of 300 arcseconds from the 
pulsar position. Beyond the distance of 300 arcseconds the emis- 
sion seems to bend sideways significantly. Since the model (see 
Sec. H]i assumes a radial symmetry, we did not extract spectra for 
any regions beyond this distance. 

In order to avoid systematic effects from the CCD borders, 
we extracted the spectra for each detector CCD separately. The 
background for each CCD was estimated using infield back- 
ground regions from the same chip. This minimizes systematic 
uncertainties on the flux, compared to a background estimation 
from blank-sky observations. The effective areas and energy re- 
sponses for each detector CCD were calculated by weighting the 
contribution from each pixel with the flux using a detector map 
from the 0.2-7 .OkeV energy band. The spectra from the different 
CCDs for each extraction region were then merged by weighting 
the responses (rmf) and auxiliary responses (arf) with the size of 
the respective region. 

For each ring, the spectra were fitted in parallel with an 
absorbed power-law model. In the individual observations that 
were used, the bad columns and the CCD borders have differ- 
ent orientations. This leads to a difference in the norm and thus 
in the flux of the fitted spectra. In the case where bad columns 
obscure a large fraction of the extraction region, a wrong flux is 
obtained. We added a selection criterion on the observations used 
for the spectral analysis of each ring to circumvent this problem. 
Only those observations were used in which the bad columns do 
not obscure parts of the particular extraction region. In this way 
we obtain the correct integral flux for each region. The result- 
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Table 3. Results obtained by fitting a power-law spectrum to the XMM-Newton data of MSH 15-52. 



Ring 

No. 



Flux 
(10~ 12 erg crrT 2 s ) 



(10 



Surface Brightness ^ 2 /dof 

17 erg cirT 2 s~'arcsecT 2 ) 



1 


1.66 + 0.02 


12.0 ± 0.3 


162 ± 4.5 


147/161 


2 


1.78 ±0.03 


10.0 ±0.3 


92 ± 3.0 


161/147 


3 


1.88 ±0.02 


2.9 ± 0.07 


9.2 ±0.2 


431/459 


4 


1.96 ±0.02 


3.1 ±0.09 


6.1 ±0.2 


356/536 


5 


2.07 ± 0.05 


2.4 ±0.1 


3.4 ±0.2 


97/156 


6 


2.24 ± 0.28 


0.4 ± 0.2 


0.5 ± 0.3 


62/107 



XMM-Newton: MSH 15-52 Spectrum - Ring 1 




Channel Energy (keV) 



Fig. 2. Spectrum and power-law fit for the first extraction re- 
gion (ring 1). Two observations were used for the analysis of 
this region (cf. Table|2]). The resulting fit parameters are listed in 
Table S 



ing number of observations for each extraction region is given in 
Tabled 

The parameters of the spectral power-law fit — absorption 
column density, photon index and normalization — were linked 
for the fitting of the different observations. The absorption col- 
umn density does not vary significantly between the different 
regions and was fixed to the value of the first ring. The fit range 
was from 0.5-9. OkeV for rings 1 to 5, while for ring 6 we chose 
a narrower fit range of 3. 0-9. OkeV, due to the insufficient statis- 
tics at low energies. As an example, Fig. [2] shows the fit of an 
absorbed power law to the observed emission of ring 1 . 

The resulting parameters for the spectral analysis of each 
ring are listed in Table [3] The data are fitted very well by pure 
power laws. Due to the increasing size of the extraction region 
and the constant area of the infield background, the statistical 
uncertainty of the spectra of the outer regions is greater, result- 
ing in a lower value of ,y 2 /dof. We obtain an absorption density 
of N H - (1.15 ± 0.03 ) xlO 21 cirT 2 using the abundance tables 
of IWilms et al.l (l2000h . This is in good agreement with results 
of previous X-ray analyses of this source (see e.g. G02). The 
extended emission of MSH 15-52 shows variation in the spec- 
tral properties with increasing radial distance to the pulsar. The 
photon index F of the fitted absorbed power law increases from 
the inner regions to the outer regions by roughly 0.5, reflecting 
a softening of the spectrum. Figure[3]shows this variation of the 
photon index. Table [3] shows the unabsorbed flux and the sur- 
face brightness in the energy range from 0.5-9.0 keV for each 
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Distance to Pulsar [arcsec] 

Fig. 3. Spectral index of the power-law fit to the XMM-Newton 
data of the different regions versus the mean distance of the ex- 
traction regions to the pulsar. A steepening of the spectrum is 
observed with increasing distance. 



extraction region as well as the x 2 an d the number of degrees 
of freedom of the fit. The surface brightness decreases with in- 
creasing distance of the extraction region to the pulsar, as can be 
seen in Fig. [4] 

For the parameter optimization of the model (see Sec.[5]l, we 
divided the spectrum for each of the rings in six intervals in en- 
ergy and calculated the flux in these bins. Only four bins were 
used for the analysis of the ring 6 region, due to the narrower en- 
ergy range. The approach of using individual flux points rather 
than only using the index and the total flux of a power-law fit 
is better suited for the parameter optimization. In this way we 
compare the data and the model results directly, without assum- 
ing power-law distributions for fits of the spectra. 



4. The Model 

4.1. Injection Spectrum 

At the termination shock of the pul sar wind, particles are ac cel- 
erated and injected into the PWN dGaensler & Slandl2006l) . In 
our model, we assume a certain particle injection spectrum at the 
termination shock and propagate the particles radially outward. 
The shape of the spectrum is chosen based on observational con- 
straints rather than on a detailed physical model, which would be 
beyond the scope of our model. We assume that the particle in- 
jection spectrum Q(E e , t) at the shock radius Rs follows a power- 
law distribution. Earlier papers on the modeling of PWN con- 
cluded that the injection spectrum should follow a broken power 
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Fig. 4. Surface brightness of the X-ray emission of the PWN in 
the 0.5-9. OkeV band as a function of the distance to the pulsar. 

law to account for the observed spectral properties at radio wave- 
lengJhsJKennel & Coroniti 1984b; Re ynolds & Cheva lier 1984; 
IVenter & de Jagerll2007l) : 



Q(E e ) = 



Go bf 



Ge 



for E e < E/,, 



for E e > E b , 



(1) 



where Qo is a normalization constant, E e the lepton energy and 
Eb the break energy of the spectrum. The parameters T\ and Y2 
are the spectral indices of the two power-law components. The 
normalization of the particle spectrum can be calculated by not- 
ing that 



f 

Je, 



Q(E e )E e dE e = rjE. 



(2) 



In this equation the parameter rj denotes the conversion effi- 
ciency of spin-down luminosity E of the pulsar into the energy 
of the leptons. The lower integration bound, corresponding to 
the minimum energy of the leptons, can be constrained by the 
observed radio spectrum and by the fact, that it has to be greater 
than the electron rest mass, E ejn [ n > 8.2 x 10~ 7 erg. The maxi- 
mum energy of the electrons is constrained by two conditions as 
specified in more detail in de Jager & Diannati-Ata'i (2009). On 
the one hand, the Larmor radii Ri of charged particles confined 
within the shock radius must be less than Rs . The parameter e is 
defined as the ratio of these two radii, e = Rl/Rs , and should be 
less than 0.5. This limits the maximum lepton energy to 



crE 



< eeK- 



c(l + ex) 



(3) 



In this expression, k is the compression ratio at the shock and 
cr is the magnetization parameter (ratio of the magnetic energy 
flux to the particle energy flux). The second condition comes 
into effect for strong magnetic fields. In this case the maximum 
lepton energy is limited by the synchrotron emission that the 
lepton radiates. Thus, the maximum energy is limited to 



<43.84fi~ 1/2 erg 



(4) 



where B is the magnetic field strength in Gauss. The lower E e >n 
obtained with the two limits is then used for Eq. [2] 



Our model introduced in this article is aimed at modeling the 
nonthermal emission of the PWN close to the pulsar. In this re- 
gion, only a young population of leptons is expected. It is there- 
fore reasonable to neglect the time-dependence of the injection 
spectrum for our model. Furthermore, we focus on the modeling 
of the X-ray spectrum and the inverse Compton (IC) spectrum in 
the very high energy (VHE) y-ray regime. For this purpose, only 
the power-law component beyond the break energy is relevant 
and we will thus assume a pure power law with E em i n - Eb in 
the remainder of the paper. This implies that r\ is only a lower 
limit on the conversion efficiency, since we do not take the lep- 
tons with energies below Eb into account in our model. 

4.2. Propagation of the Leptons in the PWN 

In the downstream flow beyond the termination shock the par- 
ticles lose energy, causing the shape of the lepton spectrum to 
change. In our model we use the approximation of a spherical 
symmetry. This implies that the propagation of the particles in 
the nebular flow of the PWN is strictly radial. Looking at the 
X-ray sky map of MSH 15-52 (see Fig. [TJ, a deviation from a 
spherical symmetry is apparent. However, compared to current 
models in which no spatial variation of the properties of the lep- 
ton population is taken into account, the radial approach should 
provide a good first estimate of the effects of spatial variation in 
the PWN. 

Assuming that the ma gnetic field in the PWN is toroidal, 
iKennel & Coronitil d!984af) showed that the product of the flow 
velocity v, the distance from the pulsar r and the magnetic field 
strength B is a constant: 



Bvr - B$vsRs - const., 



(5) 



where the subscripted parameters denote the parameters at 
the shock. This is valid for a steady-state solution in the 
ideal magnetohydrodynamic limit. The post-shocked magnetic 
field strength is then given by (IKennel & Coronitil Il984at 
ISefako & de Jagerll2003h 



Bs 




(6) 



In this equation, and for the parameter optimization, we use 
£ = k Vcr/(1 + cr) as a combined parameter for the shock com- 
pression k and the magnetization cr. 

Since we assume the approximation of a radial outflow of the 
leptons in the nebular flow, we can describe the bulk velocity of 
the leptons by a profile depending only on v$ and the distance to 
the pulsar, r. We adopt a profile of the form 



v(r)=v 5 (-^-J , 



(7) 



with an exponential index a. By using this relation and Eq. [5] 
we are able to write the magnetic field as a function of the dis- 
tance to the pulsar and thus it is possible to calculate the flow 
parameters B and v for every location in the PWN. With this 
information, we can determine the time spent by the leptons in 
different regions of interest throughout the nebula. Note that we 
neglect three aspects in our approach: the propagation of internal 
energy, the conservation of the total energy (Eqns. 5.3 and 5.4 in 
Kennel & Coroniti 1984a) and the dynamical evolution of the 
system (assuming a steady-state solution like Kennel & Coroniti 
Il984al) . The physics of these neglected aspects is contained in 
the index of the flow velocity a in Eq. [7] In the fit to the data 
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we find that the spatial variation of the X-ray and Te V spectra is 
already adequately described by this parametrization, so that an 
introduction of additional free parameters would not necessarily 
give more information about the system. 

The next step is to look at the change of the lepton spectrum 
as the particles propagate away from the pulsar wind shock. The 
particles lose energy due to adiabatic losses and synchrotron ra- 
diation in the magnetic field o f the PWN. The total ener gy loss 
of the particles is given by (cf. lde Jager & Hardindll992h : 

dE E 

= _^ V • v ± (r) - 2.368 x 1CT 3 (B ± (r)E e ) 2 , (8) 
at 3 

where the first part represents the adiabatic losses and the second 
part the energy losses due to synchrotron radiation. 



4.3. Emission Spectra 

The parameters of the lepton population and the magnetic field 
are described for all locations in the PWN by the Eqs. [5] to [8] 
With this information we are able to calculate the spectra of the 
synchrotron and inverse Compton radiation emitted by the lep- 
tons. For the calculation of the two emission processe s, we use 
the standard equations by Blumenthal & Gould (1970). 

We are now able to calculate the total emission for the lep- 
ton population at each point in the PWN. The emission from 
the model calculations is summed up to match the regions as ob- 
served in the XMM-Newton analysis (see Sec. [3}. By comparing 
the measured synchrotron spectra with the calculated spectra we 
are able to constrain the free parameters of our model. 



5. Parameter Optimization 

The optimization of the parameters of the model, as introduced 
in Section|4] was carried out by dividing the PWN into a number 
of shells. As discussed in the previous section, the values of v and 
B are known for every location in the PWN, thus enabling us to 
calculate the amount of energy lost by a particle during its propa- 
gation from the shock up to a specified shell. This, in turn, allows 
us to calculate the lepton spectrum that enters a shell, as well as 
the change in the lepton spectrum as it propagates through the 
shell. To calculate the corresponding nonthermal spectra, the ini- 
tial lepton spectrum is used. This is sufficient, provided that the 
size of the shells is chosen small enough. The modified spectrum 
is then used as the injection spectrum for the following shell. For 
the first shell, the lepton spectrum is calculated by making use 
of Eq.Q] i.e. the lepton injection spectrum at the shock. 

Since our modeling focuses only on the X-ray synchrotron 
component, we use a single power law, with E > Et,, for the 
lepton injection spectrum. The index of the lepton spectrum is 
derived from measurements of the photon spectrum close to 
the termination shock with the Chandra X-ray telescope, which 
yielded a photon index of around 1.5 (Y09). Assuming that the 
lepton population emitting this spectrum has not undergone any 
significant cooling, the spectral index of the lepton spectrum is 
then equal to 2, which we adopt for our optimizations. G02 con- 
cluded that the spectral break between the comparatively flat ra- 
dio spectrum and the steeper synchrotron spectrum is just be- 
low the X-ray band. Thus it is reasonable to assume E e ^ m i n to be 
of order 1 erg. Since only the logarithm of E e ^m contributes to 
the normalization of the spectrum, a variation of E em { n does not 
change the spectral shape and has little effect on the total flux. 

The remaining free parameters of the model are then the 
shock radius Rs , the bulk velocity of the leptons at the shock 



Table 4. Model parameters for the different scenarios as dis- 
cussed in the text. 



Parameter 


Scenario 


I 


11 


III 


Rs [pc] 


<0.1 


0.5 


0.225 


y s [c] 


1/3 


a 


no constraint 


1 


0-1 


t 


no constraint 



vs , the index of the velocity profile a, the conversion efficiency 
of spin-down luminosity into lepton energy 77 and the parame- 
ter £, which links the magnetization <x and the compression k of 
the shock. The parameters Rs and v$ are restricted by the obser- 
vations of G02 and Y09, whereas the other parameters are left 
unconstrained for the optimization. To determine the optimum 
fit to the synchrotron spectra for the different scenarios, we min- 
imized the x 1 test statistic for the XMM-Newton data points and 
the calculated model flux points. 

For PSR B1509-58, G02 and Y09 both estimated the termi- 
nation shock radius using Chandra data, but do not find consis- 
tent results. G02 favor a scenario in which a feature at 0.5 pc 
corresponds to an internal structure in the termination shock. 
They also estimated the a value for several compact knots of 
emission less than 0.5 pc away from the pulsar and concluded 
that the transition to a particle-dominated wind should occur 
at a distance less than 0.1 pc from the pulsar. Y09 analyzed a 
considerably larger data set of Chandra observations and used 
image enhancement techniques to resolve the detailed structure 
of the inner region around PSR B 1509-58. According to their 
results, the termination shock is located at a distance of 9 " to 
the pulsar, which corresponds to a termination shock radius of 
Rs = 0.225 pc. For the magnetization parameter, Y09 derived a 
value oftrsi 0.01, which is about a factor of 2 greater than the 
result obtained by G02. 

Based on the two analyses, we considered three scenarios 
with different values for Rs for the optimization procedure of our 
model. In Scenario I we assume that the termination shock is un- 
resolved in the Chandra observations by G02 and Y09 and thus 
is at a distance closer than 0. 1 pc from the pulsar. For Scenario II 
we assume a distance of Rs = 0.5 pc, while for Scenario III we 
adopt the value of Rs - 0.225 pc, as stated by Y09. 

G02 and Y09 both assume a shock velocity of vs » c/3. 
The same value was also found f or other PWN, e.g. the Crab 
Nebula dKennel & Coronitilfl9 84a). Furthermore, this shock ve- 
locity is in good agreement with the lower limit that G02 de- 
termined from the energetics of the PWN. Therefore, we also 
adopted this value for our calculations. We do not constrain the 
remaining parameters a, rj and £, except for the straight-forward 
assumption that the conversion efficiency 77 should be in the in- 
terval < 77 < 1. Table |4]gives an overview of the three scenar- 
ios and the constraints on the parameters, which we used for the 
modeling described in the next section. 



6. Results of the Modeling 

Using the model and the parameter constraints described in the 
previous section, we calculated the optimum parameters for the 
three scenarios defined in Section|5] The results of the optimiza- 
tion show that the Scenarios II and III (larger shock radius) yield 
a better fit to the data than Scenario I (small shock radius). For 
Scenario I we do not have a fixed value for Rs . We thus leave it as 
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Fig. 5. Data and model spectra for the ring 1 and 6 regions for 
Scenario II. The shaded regions mark the error bands of the 
XMM-Newton measurement of ring 1 and 6. For ring 1 the error 
band is very narrow and thus hardly visible. 
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Fig. 6. Spectral indices of power-law fits to the XMM-Newton 
data of the different regions in the energy range of 0.5-9.0keV. 
The model points shown are calculated for the best fit of 
Scenario II and use the same energy range for the fit. 



a free parameter for the simulations. The best fit for this scenario 
is obtained for a value of R$ — 0. 1 pc which is the upper limit for 
Rs (see Table 0). However, even the best fit of Scenario I does 
not give a satisfactory result. Due to the small shock radius, the 
synchrotron cooling is very efficient and leads to a strong cutoff 
in the X-ray spectra for the outer rings, which is not consistent 
with the observed emission. 

The results obtained for Scenario II yield a better fit to the 
data. In this scenario the shock radius is considerably larger. 
Thus, the leptons accelerated at the termination shock suffer 
from less synchrotron cooling up to ring 1 of the XMM-Newton 
measurement. Figure [5] shows as an example the data and model 
synchrotron spectra for rings 1 and 6 for the optimum set of pa- 
rameters found in Scenario II. Even the best-fit model spectrum 
deviates significantly from the XMM-Newton spectrum, but this 
is expected due to the simplifications of our model. The general 
trend, however, reproduces the spectral shape and the observed 
flux level. The spectral index with increasing distance of the ex- 
traction region from the pulsar is plotted in Fig. [6] The best-fit 
model is able to reproduce the variation of the spectral index. 
As shown in Fig. [7] the change in surface brightness with in- 
creasing distance of the extraction regions to the pulsar is also 
reproduced. For Scenario III we also get a good fit to the data. 
Figure [8] shows the spectra for the optimum set of parameters 
found in Scenario III, again for rings 1 and 6. The spectral in- 
dex and the surface brightness with increasing distance can be 
seen in Figs.|9]and[T0] The shape is also reproduced well for this 
scenario. 

The results obtained for Scenarios II and III constrain the 
range of the model parameters. Since the parameters are corre- 
lated, it is only possible to state parameter ranges. For Scenario II 
we found that the index of the velocity profile is in the range of 
a = 0.4 - 0.6. The conversion efficiency is greater than 0.3 for 
this scenario and £ ranges between 0.3 and 1.3. The compres- 
sion ratio k varies between 1 and 3 for relativistic shocks. We 
can thus translate the range of £ to the result that <x > 0.01, 
which is a factor of 2 more than the lower limit derived by G02. 
For Scenario III, we found the same range for a and 77, but a nar- 
rower lower limit on the magnetization parameter, <x > 0.005. 
This is lower than the limit that Y09 found in their estimate 
based on the equipartition assumption. In summary, we can state 
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Fig. 7. Variation of the surface brightness for the XMM-Newton 
data and the model in the energy range of 0.5-9.0 ke V. The model 
points shown are calculated for the best fit of Scenario II. 



that our results support a shock radius of the order of 0.2 - 0.5 pc . 
However, we are not able to favor one scenario above the other. 
For the conversion efficiency, we found a lower limit of 77 > 0.3 
for Scenario II and III. This agrees ve ry well with t he res ult of 
the time-dependent one-zone model by Zhang et al. (2008). The 
magnetic field esti mates for MSH 15 -52 range between 8//G 
(G02) and 25//G dZhang et al.ll2008l) . The spatial evolution of 
the magnetic field with distance for our model can be seen in 
Fig- EH Our predictions for the magnetic field in the PWN agree 
with the estimates from the other authors. The lower limit of 
<x > 0.005 on the magnetization parameter agrees very well with 
the observational results obtained with the Chandra satellite by 
G02 and Y09, which both state a lower limit of the same value 
as we derive with our model. 



7. Implications for the TeV y-ray Emission 

Based on the optimized parameters described in the previous 
section, we calculated the IC emission for our model to com- 
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Fig. 8. Data and model spectra for the ring 1 and 6 regions for 
Scenario III. The shaded regions mark the error bands of the 
XMM-Newton measurement of ring 1 and 6. For ring 1 the error 
band is very narrow and thus hardly visible. 
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Fig. 10. Variation of the surface brightness for the XMM- 
Newton data and the model in the energy range of 0.5-9. OkeV. 
The model points shown are calculated for the best fit of 
Scenario III. 
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Fig. 9. Spectral indices of power-law fits to the XMM-Newton 
data of the different regions in the energy range of 0.5-9.0 keV. 
The model points shown are calculated for the best fit of 
Scenario III and use the same energy range for the fit. 



pare it to the observational TeV y-ray data. MSH 15-52 is an ex- 
tended sour ce in VHE y-rays and w as detected by the H.E.S.S. 
experiment (lAharonian e t al. 2005). The emission in the VHE 
energy band extends beyond the central part of the PWN that is 
considered in our model. Thus, we used the published VHE ex- 
cess events map to rescale the spectrum to the central 300 arcsec 
region around the pulsar. The excess map, smoothed with a 
Gaussian with the size of the point spread function (a;0.07°, is 
shown in Fig. Q~2] The ratio of excess events between the two 
regions used in the H.E.S.S. publication and for our modeling 
is 3.42. The rescaled normalization of the y-ray spectrum at 
1 TeV is then (1.67 ± 0.06 stat + 0.41 sys )x 10~ 12 TeV" 1 ctrT 2 s _I , 
where we rescale the errors with the same factor. Our rescal- 
ing assumes that there is no significant spectral variation for the 
H.E.S.S. spectrum, i.e. the photon index of the measured index is 
equal to 2.27 ± 0.03 sta t ± 0.2 sys throughout the extraction region. 
Another important aspect to consider is the size of the H.E.S.S. 
PSF Since the extension of its 68% containment radius is of 
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Fig. 11. Spatial evolution of the magnetic field with distance 
from the pulsar. The two solid lines denote the upper and lower 
range of values of the magnetic field strength for Scenario II 
(based on the range for a and ^ that is found in the optimiza- 
tion). The two dashed lines mark the upper and lower limit of B 
for Scenario III. 



the same order as our extraction region (/?68~0.1deg), we ex- 
pect some contamination in our extraction region from the outer 
regions of the PWN. About 30% of the emission from inside 
the PSF will fall outside of the region defined by the rings 1 to 
6. However, a similar amount from the emission outside the PSF 
region will fall into the central PSF region (i.e. ring 1 to 6). These 
effects will cancel to a first order, although a systematic steep- 
ening in the spectrum is expected if the TeV photon index also 
steepens with increasing radius, as is predicted by our model. 
Figure [13] shows the spectral index for the IC emission with in- 
creasing distance of the region from the pulsar, where it can be 
seen that the index steepens by Ar w 0.15. 

Figure [14] shows a plot of the spectral energy distribution 
(SED) of MSH 15-52 in the X-ray and TeV y-ray band. The 
presented experimental data in the X-ray band are the sum of 
the emission of the shells from the XMM-Newton analysis (see 
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Fig. 12. Smoothed H.E.S.S. excess map of the VHE y-rays from 
MSH 15-52 ( Aharonian et alj2005l) . The map was smoothed by 
a Gaussian with the size of the point spread function (a:0.07 o ) 
to reduce statistical fluctuations. The extraction region for the 
published H.E.S.S. spectrum is marked with a dashed black line, 
the white solid line denotes the 300 arcsec region used for the 
modeling in this paper. The position of the pulsar is indicated by 
a white cross. 



Sec. |3). For the TeV range, the rescaled H.E.S.S. spectrum is 
shown. The model SED was obtained with the best fit param- 
eters of Scenario II (Rs- 0.5 pc) for a conversion efficiency 
of T] = 0.4. For this set of parameters, the fitted spectral in- 
dex to the IC emission in the 0.3-30 TeV energy range is 1.8 
and thus slightly below the statistical and systematic errors of 
the measured H.E.S.S. spectrum. The normalization of the fit is 
3.1x 10~ 12 TeV ' crrT 2 s and hence also below the lcr er- 
rors of the measurement. The ratio of the predicted flux above 
0.3 TeV for this set of parameters to the observed flux is 0.3. 

The change in surface brightness of the IC emission is shown 
in Fig. [15] The flux was calculated in the energy interval of 0.3- 
30 TeV and for the same ring regions as used in the X-ray anal- 
ysis. Unfortunately, with the angular resolution of the H.E.S.S. 
experiment it is not possible to test the model against the data. 
Future Imaging Air Ch erenkov Telescopes like the Cherenkov 
Telescope Array (CTA, [Wagner et al.l l2009) will have a resolu- 
tion around one arcminute and will allow for a test of these pre- 
dictions. Compared to the X-ray surface brightness, which drops 
by three orders of magnitude from ring 1 to ring 6, the drop in 
the surface brightness of the IC emission is considerably smaller. 
However, from looking at the VHE sky map (Fig. [T2l one can 
see, that the emission must have a large component which is not 
described by our model. This might, for example, be due to an 
older population of electrons. 



8. Conclusion 

We analyzed XMM-Newton data of the SNR MSH 15-52 to de- 
rive the spatially resolved spectral parameters of the inner re- 
gion of the source. For this analysis we extracted spectra from 
six annuli centered on PSRB1509-58. A steepening of the X- 
ray spectrum with increasing distance from the pulsar is ob- 
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Fig. 13. Spectral index of a power-law fit to the IC emission cal- 
culated in the model with increasing distance from the pulsar. 
The fit range was from 300 GeV to 30 TeV. Plotted is the result 
of the model with the best fit parameters for Scenario II for a 
conversion efficiency of r\ — 0.4. 
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Fig. 14. Spectral energy distribution of MSH 15-52 from the ex- 
traction regions as marked in Fig. [T2] The XMM-Newton and 
the H.E.S.S. data are marked as solid lines, the model results for 
the best fit values of Scenario II with a conversion efficiency of 
rj - 0.4 are drawn as dashed lines. 



served. The surface brightness in the XMM-Newton range drops 
by three orders of magnitude from the inner region close to 
PSRB1509-58 up to the last ring used in our analysis (at a dis- 
tance of 300 arcsec from the pulsar). The spectra of all rings are 
fitted well with power laws. 

We then introduced a numerical model to describe the spa- 
tial evolution of the lepton population and the magnetic field 
within the PWN. Our model is based on the parameters of the 
termination shock. We fitted the calculated model spectra to the 
XMM-Newton data to constrain the parameters of our model. 
The results of our optimizations suggest a termination shock in 
the range of 0.225 to 0.5 pc (Scenarios II and III in this work), 
which are the values derived from different analyses by G02 and 
Y09. However, the results of the optimizations of our model do 
not favor one scenario over the other. For the magnetization pa- 
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Fig. 15. Surface brightness of the IC emission with increasing 
distance from the pulsar. The surface brightness was calculated 
for a power-law fit to the data in an energy range of 300 GeV to 
30 TeV. Plotted is the emission from the model calculations with 
the best fit parameters for Scenario II for a conversion efficiency 
of j] = 0.4. 



rameter, we arrive at a lower limit of 0.005, while our lower limit 
on the conversion efficiency is 0.3. 

Using the model parameters obtained from the fit to the X- 
ray data, we calculated the IC emission from the lepton popu- 
lation in the PWN. We rescaled the published H.E.S.S. data to 
compare the total IC emission from our model with the observed 
emission. The results show that our model reproduces the right 
order of magnitude for the observed IC flux, but differs slightly 
in the spectral shape of the emission. However, we did not take 
into account contamination effects for the VHE photons due to 
the PSF. A steeper photon index outside the region used for our 
model would lead to a systematically steeper spectrum inside the 
region due to PSF spillover effects. This effect was not taken into 
account in Fig. [14] Additionally, we also calculated the spectral 
index and the surface brightness with increasing distance of the 
regions from the pulsar, which may be tested with future VHE 
y-ray instruments like the CTA observatory. 

The total energy that is emitted by the leptons in the in- 
ner region of the PWN via synchrotron and IC emission in our 
model can be obtained by integrating over the SED (see Fig.fPfl). 
The integrated emission of 7.6 x 10 35 erg s -1 amounts to roughly 
4% of the current spin-down luminosity of PSRB1509-58. 
However, when making this comparison one has to be aware 
that the integrated emission comes from a population of leptons 
with different lifetimes and cannot be directly linked to the cur- 
rent spin-down luminosity. Furthermore, there is also significant 
emission at X-ray and TeV energies outside the central 0.1° re- 
gion considered in our modeling approach, which increases the 
percentage of lepton energy converted to radiation (see Figs. Q] 

andO- 

In summary, we show that our radially symmetric model of 
the PWN MSH 15-52 readily describes the nonthermal emission 
seen in the X-ray band. Despite the simplifying assumptions of 
our model, the change in spectral index and surface brightness 
can be reproduced. The model is also able to predict the flux 
level of the observed emission in VHE y-rays, however, the ex- 
act shape of the spectrum is not reproduced. For the future, an 
application of the model introduced in this paper to other PWN 



should be interesting to obtain a broader view on the underlying 
physics in PWN. 
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